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Abstract 

We have analyzed the dependences of the threshold value of amplitude of 
linear density fluctuation collapsed at the current epoch, 6c, and its overdensity 
after virialization, A^, on matter density content, 3D curvature parameter 
and cosmological constant. It was shown that in models with negative or 
zero curvature {k < 0) and positive cosmological constant (A > 0) despite 
smaller rate of perturbation growth, 6c appears to be smaller than its canonical 
value 1.686 [k = A = 0) due to the longer duration of amplitude increase. 
However, in models with positive 3D curvature, 6c can be larger when the 
matter density parameter is small. In this case critical amplitude 6c approaches 
infinity when converges to its small critical value. Though the range for the 
threshold values of perturbation amplitude is quite narrow - 1.55 < 6c < 1.75 
in the region of parameters —0.4 < ilfc < 0.4, < Qa < 1, 0.1 < Qm < 1, 
the difference in the concentrations of rich clusters of galaxies (as calculated 
within the Press-Schechter formalism framework) with real and with canonical 
values of 6c reaches, and for some models exceeds, 30%. The range of changes 
for the overdensity after virialization, Ac, is considerably wider for the same 
region of parameters: 60 < Ac < 180. It results into difference up to ~ 
40% between the X-ray temperatures of gas, as calculated for these values 
and for the canonical value of Ac = 178. Also we have found analytical 
approximations of the dependences 6c{^m,^A) and Ac{i^m.,^A) with their 
accuracies for above mentioned region of parameters being no worse than 
0.2% and 1.7% respectively. 
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1 INTRODUCTION 

The mass function and X-ray temperature function of rich clusters of galaxies are 
important characteristics of the large-scale structure of the Universe. The studies of 
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rich clusters are so intensive due to their eminent properties: a) they are the largest 
gravitationally bound systems in the Universe; b) the observational selection crite- 
rion for these structures is quite evident enabling one to draw the catalogues of these 
objects with the highest completeness, at least within the radius of ~ 300/;,~^Mpc, 
within cones ~ 30° around the poles of our galaxy; c) space and velocity distributions 
of galaxies in the clusters demonstrate their dynamical equilibrium - hence indicate 
to the possibility the main theorems of dynamics of gravitational systems to be ap- 
plied for determination of their characteristics at current state; d) Press-Schechter 
formalism and its modifications [1-9] connect these physical parameters with fea- 
tures of the power spectrum of initial density perturbations. This power spectrum 
is sensible both to dynamical parameters of the Universe (Hubble constant, 3D-space 
curvature, cosmological constant), which determine the rate and acceleration of its 
expansion, and to the density parameters of different matter components (baryons, 
photons, neutrinos and hypothetical particles which compose the hidden mass of 
the Universe) - that enables one to determine these parameters on the basis of the 
observed data on the large-scale structure of the Universe. 

The latest experimental programs dealing with detection of the microwave back- 
ground temperature fluctuations on different angular scales (COBE, Boomerang, 
MAXIMA, DASI etc.) and more complete catalogues of galaxies and galaxy clusters 
(Sloan Digital Sky Survey, 2dF Galaxy Redshift Survey etc.) allow us to determine 
main cosmological parameters and accuracy of such determination. The enhanced 
experimental accuracy requires more precise theoretical description of the observed 
phenomena. Press-Schechter formalism describes dependence of the mass function 
of rich clusters of galaxies on the mean matter density of the Universe (p), disper- 
sion of density fluctuations on the scale of clusters {<Jm) and threshold value of the 
amplitude of density fluctuation collapsed at the current epoch {6c): nps{M) = 
F{M, p, aM)6c/cTM exp{-5^/2alj), where F = -y/^p/M^ \d Inau/d lnM\. Ecke 
et al. (1996) |H had shown that the observed data of those time were well ap- 
proximated by such dependence with 5c = 1.686, as could be inferred from the 
dynamics of the collapse of spherically symmetrical cloud for the flat model of the 
Universe without cosmological constant. The conclusion drawn by these authors 
regarding the weak dependence of the 6c value on 3D-space curvature and cosmo- 
logical constant paved the way for the canonical value 6c = 1.686 to be used at 
arbitrary values of 3D-space curvature and cosmological constant without the anal- 
ysis of errors considered. However, the data of the apparent magnitude - redshift 
testing for Supernovae la become available since 1998 (see references in [10-12]) and 
suggest that the Universe is expanding with acceleration. Thus the most plausible 
model of the Universe rather strongly deviates from the canonical model with a zero 
3D curvature and zero cosmological constant. Most probably, the matter density 
(dark matter plus baryons) is only 20-30% of the critical density value, and the rest 
70-80% is a dark energy (cosmological constant, vacuum energy, or quintessence). 
This requires more careful approach to applying of the Press-Schechter formalism 
for determination of concentration, mass function and X-ray temperature function 
of rich clusters of galaxies in different cosmological models for their comparison with 
observed data. 

The paper is aimed to: 
- analyze the dynamics of the collapse of spherical dust-like cloud in models with 
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arbitrary values of 3D curvature and cosmological constant; 

- to study dependence of 5c, the threshold value of the initial density perturbation 
amplitude, on the above mentioned parameters, as well as dependence of the final 
density fluctuation after establishment of dynamical equilibrium, Ac; 

- construct approximations of these dependences on the values of the 3D-space 
curvature and cosmological constant; 

- study deviation of the mass function and X-ray temperature of rich clusters of 
galaxies employing real and canonical values of 5c and Ac (1.686 and 178 respec- 
tively) for the models being analyzed in modern literature. 

The analyses of the dynamics of the collapse of spherical dust-like cloud for 
models with zero 3D-space curvature and arbitrary values of cosmological constant 
as well as for generalized models, with non-zero curvature, were carried out in [4lfT3]. 

In this paper we study the dynamics of the collapse of spherical perturbations for 
models with arbitrary curvature and cosmological constant parameters and analyze 
the influence of the dependence of critical amplitude on these parameters for the 
calculation of mass function of rich clusters of galaxies. Aiming at this, we will find 
solutions of appropriate equations, which describe the development of density per- 
turbations at all evolutionary stages from initial linear one up to nonlinear collapse 
and virialization. 



2 The dynamics equation for a dust-like cloud when 

A 7^ and /c 7^ 

In order to analyze the development of spherically symmetrical perturbations in the 
Universe with a non-zero cosmological constant (A 7^ 0), we consider the dynamics 
of dust like cloud in the comoving synchronous gauge. Note that similar problem 
was first solved by Tolman in 1934 for the case A = We will consider 

the dynamics of spherical cloud before the appearance of first counterflows, that 
is trajectory intersections of single particles (spherical shells). At this stage, the 
dynamics of matter is well described by the continuum approximation, that is by 
the energy-momentum tensor of perfect fluid 

Tik = {e + P)uiUk - gikP, («, /c = 0, 1, 2, 3), 

where e = pc^ is energy density, P is hydrodynamic pressure, Ui - 4D velocity 
components. 

The continuum approximation with P = allows to consider perturbations in 
the gauge, which is simultaneously synchronous and comoving {ua = 0, a = 1, 2, 3). 
In spatial spherical coordinates R, 9 and ip the metric of such Universe would take 
the following form [1^1 

ds^ = dr^ - e^(^'")rfi?2 _ r\R,T){de^ + sin^ Odip^), (1) 

where r = ct is the time in units of length. Einstein's equations with cosmological 
constant are so 

^ik — -Qik^ — ^J-ik + Qik^^- 
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Calculating the tensor components Rik for the metric JH) and substituting them into 
Einstein's equations, and assuming T^p = for our case (where a and /5 = 1,2, 3), 
we obtain the following set of equations: 

g-^ 1 
-(2r"r + - rr'X') + —(rrX + + 1) = SnGp + A, (2) 

2rr + P + l- r'^e"^ = Ar^ (3) 

A + 2- + ^ + — - — (2r" - AV) = 2A, (4) 
r 2 r r 

Xr' - 2f' = 0. (5) 

All other components are identically zeros. By integrating the latest equation, the 
following relation is obtained: 

i-fiRV 

where the value f{R) satisfies the condition 1 — f{R) > 0. Substituting it into (jnj 
results in 

2rr + - Ar^ + / = 0. (6) 
Integrating of this equation gives 

- -^-^ - -r' = -/, 7 
r 3 

where F(R) - time independent value. This equation can be interpreted as an 
energy conservation equation for unit mass particle moving in the potential field 
U = — (A/6)r^ — F/(2r) with total energy equal to — //2. Multiplying ^ by r'P 
and using ^ and (0), we obtain 



SvrGp 



F' 



Carrying out the integration of the latest equation it becomes 

R 

F{R) = 8ttG J pr'^r'dR = 2Gm, (8) 



where m - total mass of the sphere with comoving radius R. Since integrals are 
taken over the hypersurface of constant time, so dr = r'dR. By substitution p = 
p(l + 5(i?,r)) (where r) = {p{R, t) —p(t)) /p{t) andp(r)-mean matter density 
in the Universe) into equation (jHI) we obtain 

F{R) = ^pr%l + 6r) = Hln^ Q)' (1 + 5,), (9) 

where 

3 ^ 



3 f 

6riR,T) = S{y,T)y^dy. 
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Also we have used here the definition for the mean matter density in the Universe 



where is the current value of the Hubble constant, - the density of dust- 
like matter in units of the critical density of the Universe, air) - scale factor of 
cosmological model with metric 



ds = dr —a ir) 



dR'^ + x\R){d9^ + sin^ed^^' 



Therefore 



or 



Wr) = ^f-V-1, (10) 



r,' \ ^ /a 



6r{R,T) = {6r{R,n) + l)i-) - -1, (11) 



where the index i denotes fixed values at a given period of time Tj. In the absence 
of perturbations, the components of metric take the form 

riR,T)=aiT)xiR). 
In this case, the set of equations ((SHSI) will be reduced to the following 

a' - - 3 « = -k, (12) 

2da + - Aa^ + k = 0, 

x"{R)x{R)-x'\R) + i = o, 

where k = {1 — x'"^ [R)) / {R) - curvature parameter. The real solutions of the last 
equation are 

R) / oj 
R 

sh{ujR)/uj 

where u = h/Rq, and Rq - radius of the Universe. With them the curvature 
parameter takes values 




For closed models (k = 1) Rq is real number, in the case of open models (k = —1) 
radius of the Universe is imaginary number. The set of equations (0), (fTI|) and (fT2|l 
allows us to describe the evolution of perturbations from early linear stage up to 
late nonlinear one. The initial conditions are following: 6r,i = 6r{ai) and rj = r(aj). 
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3 Linearization of Tolman's model equations and 
their solutions 



At first, let us consider this set of equations for the case of small perturbations, 
namely for S 1. Then the metric components (pj) may be written in the following 

way 

gA(«,r) ^ ^2^^)^^ ^ hu{R,T)),r{R,T) = a(r)x(/?)(l + h22{R,r)). (13) 
Employing (0) and one may represent the equation Q in the form: 

2fr H p(l + 6r)r^ = 0. 

3 3 

For small perturbations {hu <^ 1, /i22 ^ 1) the last equation in the linear approxi- 
mation acquires the form 

a^/i,22 + 2da/i22 H —pSr-a^ = 0. (14) 

3 

In order to find the relation between the /i22 and 6r, let's insert (fT3|l into (fTfll and 
confine ourselves by the 1-st order values 

^ F(i?) / _ _ 

mn^x'iR) [ 2 ; • 

Switching onto the case of non-perturbed metric {h22 = 0), one obtains 

6r = F{R)/{Hln^x\R)) - 1, 

with (5r = 0. The presence of density perturbations in the matter distribution 
inevitably leads to respective perturbations in the metric, that is /i22 7^ 0. Therefore, 
the amplitude 5r for /122 = must to be identically 0. Hence 

F{R)=Hln^x\R)- (15) 

In other words, the analysis of perturbation behavior is carried out by comparing 
space and velocity distributions of the same particles in the regions that contain 
equal mass m within the same comoving R. The equation for the growth of matter 
perturbations (fTUI when (fT^ is assumed, becomes the next 

This equation is often used for the correction of linear power spectrum for the 
nonlinear evolution of the perturbation [12]. Taking it into account one may obtain 

h22 = (16) 

In fact, this condition let to single out from the arbitrary initial conditions those ones 
which determine a scalar perturbation mode - of the most interest in cosmological 
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studies. The substitution of the latest equation into (fT^ provides the following 
linear equation for the evolution of the small perturbation (see also p[7]) 



6r + 2-6r - AnGp6r = 0. (17) 
a 

The general solutions can be expressed by the relative difference in background 
densities of the two similar models which differ in integration constants 

1 da 

0, 



a daj 



where aj - integration constant of the cosmological equations for background. 

To find the solutions of (fT2|l we need the initial value of the scale factor = a(rj) 
for some arbitrary time Xj. In particular, for the start time r = we fix the value 
of scale factor a(r = 0) = 0. Due to this we can obtain the value of scale factor 
from (fT^ for any time and in particular for this moment Tq. The current time Tq 
is obtained from Hubble constant. This constant is denoted as Hq at current time, 
and the scale factor as oq. In the case of non-zero curvature 



ao 



k 



Values of pg and A at present time are defined through physical coordinates and not 
dependent on value of Oq. In the case of zero curvature the scale factor is indefinite 
and could by taken deliberately. Apparently, the physical sense is contained in the 
ratio a(ri)/a(r2), where Ti and T2 any moments of time you choose, not in a base 
value of scale factor. Let us reduce the scale factor to its present value oq. So, when 
noting this rate as a = a/oo, the equation (fT^ takes the form 

^2 87rG'__2 A_2 r 
a pa a = —k, 

where k = k/a^ takes continuous values. Below we denote k by k, considering it 
as continuous quantity, and a by a, putting scale factor at the present time equal 
to unity. The quantities tq and Oq, or, what is the same, tq and k, are integrals of 
cosmological equation. Thus, we obtain two solutions of equation (fT7j) 

1 9a Ida 

a otq a ok 

According to the above definition 

piT,k,To + dro) - p{T,k,To) p{T,k + dk,To) - p{T,k,To) 

^r,l ^7—, ^ , 0r,2 ^7—, ^ • (18) 

Using the integral of the equation (fT2|l . the last expressions can be represented as 

a ' a J 
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where X = ^a^ + ^^pa? — k. We are interested in the growing solution 5^,2 so far. 
Let's introduce the designation jT2| 

D[a) = ^iS^x'/^ Jx-'/^da (19) 
^ 

Then the exact linear solution for the increasing mode of matter density perturbation 
is 

- Ir^*"'- <^»' 

In view of the fact that perturbation at the initial stages was very small, one may 
assume 

6r{a) - 6l{a). (21) 
When a-Cl, one may expand (f2fl|) into series over a small parameter a, so that 

(5,(a) ~(5^(a) ~ A(a + 0(a')), (22) 

where 

A=M^. (23) 
D{ai) 

4 Parameters Sc and Ac, and their dependences on 
A and k 

One of the important values used in cosmology for analyzing conditions of the large- 
scale structure of the Universe formation is the initial amplitude of matter pertur- 
bations which collapse at the given moment of time. It is called threshold amplitude 
and denoted by 6c and depends both on the parameters of cosmological model and 
on the time moment of collapse Tcoi {Sc = Sc{^m,^A,Tcoi))- In the framework of 
linear theory of perturbations the increasing solution is described by the expression 
(f2n|l . To find the solution that would describe the case for 5 > 1, let's use the 
fact that the increasing mode of scalar perturbation is determined by the difference 
between the models with different curvature and the same start time Sr = 5r,2 (fT8|) . 
This makes possible to extract the increasing part out of the general solution (|TTHl 
by use of the time synchronization of integrals of equations and (fT^ respectively 

n r(a) 

da f dr 



^a^+'-^pa^-k i l^r^ + L^f 



^0, / ^ = = r - To. 



Equating them gives 



a 

I 



r{a) 

da r dr 



(24) 



At the moment when perturbation turns around, i.e. when r = 0, r reaches its 
maximal value r = r^a, where rta = r{ata) and can be obtained from the equation 

A 



3 < - frta + F = 0. 
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At the same moment the right-hand integrand in the equahty (|24|1 takes an infinitely 
large value, but the integral remains an finite quantity. So, let's rewrite the to 
the form 



a.ta . rta . 

da r dr 



where PV means principal integral value. 

As we already mentioned, the equation is equivalent to the energy conser- 
vation equation for the massive particle motion in given spherically symmetrical 
potential field. It is well known, that in the case of negative total energy (that is 
when / > 0) such motion is symmetrical in time with respect to f = moment, 
thus 

where Tcoi - moment of the collapse of the cloud. We have assumed here that the 
expansion of cloud started at zero moment of time. This equality can be rewritten 
in the following way: 

^col , f'ta , 

da f dr 

2PV 



where acoi = aijcoi)- Taking into account (fTK| . one obtains 

where = A/{3Hl), = 8nGpJ{3H^), fi^ = -k/Hl = -f/{Hlx\R)), 
X = r/x{R), xta = rta/xiR)- 

Setting the time of collapse and using the last equality, one can find the value 
Qf for the perturbation collapsing at the moment acoi- Substituting (fT3ll with (|T6| 
and (fTH| into (jlj) and confining ourselves to the 1-st order quantities: 

aa-+[^a + ^pa ' ) = - - k ■ 



Considering the equality only at small a, substituting (f^ into this equation and 
writing out components at zero power of a gives: 

Hence 



A 



Taking into account (fTTl . ^K)^ and (pTHl we obtain 



6Ua) = ' : '^ D{a). (26) 
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Thus, the initial amplitude of increasing mode of the perturbation is determined 
by the difference of curvature parameters of the background 3D-curvature parameter 
flk and local perturbation Qf one. Fixing the time of collapse Tcoi in arbitrary way 
and determining the value Qf from (f25|l . we substitute it in (f26|l and obtain the 
sought value of 6c 

Sc{acoi,a) = —-^ D{a), (27) 

where Qcoi - scale factor at the moment of collapse, a - scale factor at the moment for 
which recalculation of the critical amplitude is done according to the linear theory 
of perturbations. In the case of flk = 0, after appropriate redesignations, the above- 
stated formula is equivalent to the similar value obtained by [ij. For the simplest 
case of flm = 1 and fifc = it coincides with the well known one from literature (see 
for example [HI ED]), since, as it may be easy shown, the critical amplitude will 
take on the form 



3 /37r\2/3 



5 V 2 / ami 



a ,_ a 



1.68647- 



(^col ^col 

and at a = a^oi it has the canonical value. 

Fig. [U shows dependence of 5c on Vt^ for the models with fixed values of cur- 
vature parameters VLk (a) and cosmological constant VL^ (d). In the models with 
positive values of cosmological constant the condition of positive curvature of local 
perturbation is not sufficient for its ever collapse since the positive cosmological 
constant counteracts to deceleration at the stage of expansion of perturbation. The 
local curvature (//x^(-R)) needed for the collapse must exceed some critical value 
determined from the condition that there exists the moment of turn around of per- 
turbation region where f = 0. 

The magnitude P according to ([Zj) acquires the minimum at r = Vcxtr = (^m/ (SI^a))^/^. 
Then 



•2 -2/ 

min ' ' 



extr 



In order the turn around moment to exist there must be r^j„ < (only at this 
condition there is a moment when r = 0). That means 

As it was already shown above, the initial value of perturbation is determined 
by the difference between the background curvature and the local curvature of the 
perturbation region (see (f26|) ). For the models with negative curvature {Qk > 0), 
the initial value of the density perturbation amplitude of collapsing cloud always 
must exceed some critical value, because the collapse region must have a positive 
value of curvature. Moreover, the larger flk is the larger value of the amplitude 
of perturbation takes at the turn around of its expansion, as calculated within the 
framework of linear theory (see (f27jl and FiglUb and e) . Among background models 
with a positive curvature {Qk < 0), one can find models with considerable slow-down 
of their expansion. Moreover, the perturbation amplitude needed for collapse can be 
infinitesimal (FigHjD). However, after recalculation according to the linear theory for 
the present time (Figtl^), the initially small perturbation amplitude tends to large 
values. This is due to the fact that these models have long period of time when 
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the rate of expansion of the background model was close to zero d ~ (Fig. [2^). 
Both for the perturbation region and background, the minimum value for the rate of 
expansion, according to (fT^ . is reached when a = Uextr = i^m/ C^^a))^^^ ■ Similarly, 

In contrast to the perturbation region, the smallest rate of expansion for back- 
ground models with VLm < 2f2A {aextr < 1) have to be greater than zero. Otherwise, 
we would observe now a stationary or collapsing model of the Universe. At the mo- 
ment of the minimum rate of expansion, as it could be easily shown, the acceleration 
equals zero. Nearly zero rate of expansion and its invariability with time leads to 
the appearance of a long-lasting stage in the evolution of the Universe, where the 
value for the scale factor is almost constant. According to (fTTjl . perturbations at 
this stage grow exponentially 

~ e , 

where a = ^/4'KGpQa^^/r , a At - duration of a considerable slow-down in the 
background expansion (in units of the Hubble constant). The closer the minimum 
rate of expansion is to zero (see Fig. [2^), the longer stage is and the increase of 
amplitude is greater (see Fig. Wp)- The characteristic time for this stage can be 
obtained on the basis of (fT2|) 

T e, 

where e = Aa/a - an allowable relative change in the scale factor. At the infinites- 
imal value of e and nearly-zero value of dmm the time of exponential increase of 
perturbations At may take arbitrary large values. Correspondingly, an increase of 
the perturbation amplitude for this period of time is large as well. Thus, cosmolog- 
ical models with a substantial slow-down in expansion are gravitationally unstable. 

In models with flk > such slow-down is essentially less, resulting in an in- 
verse relation between Sc{acoi) and Sc{ata)'- Sc{acoi) / Sc{ata) decreases with decrease 
of flm- Such behavior is caused by decrease of the rate of perturbation growth at 
the final stage of the collapse of perturbations due to an increased influence of the 
cosmological constant and curvature of the Universe (Fig. E]). 

Dependence of the critical amplitude 6c = 5c(^coz) recalculated to the moment of 
collapse, allows to evaluate a change of influence of the cosmological parameters on 
the dynamics of the growth of perturbations in the expansion Universe (Fig. |4K). 
It is clear that the influence of the cosmological constant on Sc with increasing of z 
diminishes as it is also substantiated by graphs. 

Dependence of Sc{acoi = 1, a = 1) on flm and can be expressed by the following 
approximate formula: 

6c = A- B\n{nm + C), 

where 

A = 1.687 + 0.021^]A, 

B = 0.01568(fiA)^-^^^ - 0.03061, 

C = 0.00195 - O.OllSS^A + 0.00382(fiA)^- 

The numerical values of the coefficients are found by the method of best-fit ad- 
justment. The accuracy of approximation for different parameters of the model are 
shown in Fig. Eti- 
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Note that the real amplitude of perturbation 6, calculated according to (fTTH) . is 
infinitely large at the moment of collapse. In fact, at the final stage of collapse of a 
cloud the energy is redistributed with subsequent equilibrium due to the short-range 
particle interaction and a finite target parameter. At very small scales, the gravita- 
tional force is counteracted by the short-range interaction that is determined by the 
nature of the particles in a cloud. At large scales, the key role in formation of sta- 
tionary structures plays the balance of inertial and gravitational forces. According 
to the virial theorem, a system of particles, approaching its dynamical equilibrium, 
has the following kinetic energy per unit mass [IQJ 

i i 

In our case: U = U\ + Up = —^r"^ — ^, so 

6 2 2r ^ 2 ^ 

The moment when this equality is true we shall call the moment of virialization and 
denote with an index 'vir\ Let's make use of the conservation law expressed in the 
form: 

T+UA + Up = UA,ta + Up,ta. 

At the moment of virialization, the last equality takes on the form of [T9] 
or another form 



where x^ir = r^ir/x{R)- 

It would be reasonable to introduce here a quantity to characterize the relative 
divergence of density taking into account virialization process. This quantity, called 
as overdensity after virialization, is defined in the following way: 

Pvir / _ ^^ ^tuHq 



Ar uti / 
c — [O'col 



Pc ^vir-^'^^^col) 

In the case of the canonical model {VLa = 0, Qm = 1, ^fc = 0) 

Ac = ISvr^ ^ 178. 

Dependence of Ac on cosmological parameters is shown in Fig. El and Fig.Eb- At 
scales of stationary clusters of galaxies, the density of the energy that is connected 
with the cosmological constant is far less than the energy of the matter, hence it has 
a subtle influence on the final value of the radius r^j^- In this case, we may assume 
with a high certainty that r^r ~ rta/2. Thus, the final amplitude of the virialized 
perturbation is completely determined by scale of turn around Xta = ^ta/xiR)- 
When the relative density of the matter VLm increases, the perturbation amplitude 
grows nearly proportionally to it. Although for the models with small values of Vtm 
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this linearity is broken because of the increase in the role of cosmological constant 
(see Fig. and Fig. EJd). 

The behavior of overdensity Ac could be explained by analyzing the influence 
of every parameter on its dynamics. As far the matter component f2m makes the 
process of expansion more slow, so the value of should be decreasing when this 
component's content increase (one of the components, either or ^Ik we will assume 
to be fixed). So that the value Ac{acoi = 1) = ^m/a^^ir '^iU be growing. In the same 
way, since component accelerates the expansion, will be growing. This leads 
to decrease in amplitude Ac when Qm or is held fixed. 

As it was shown above, the relative amplitude of perturbation depends linearly on 
the difference of curvatures, 1)^ — fi/. The local curvature parameter, depends on 
parameters of background model for the fixed collapse moment. However, because 
of slow variation in Qf against flk at the flm or fl\ being fixed, the change in value 
of — f2/ will be determined by the change in value of curvature parameter Qk- 
Therefore, if flk increase, then Ac do the same. These relations of Ac dependence 
to the cosmological parameters are illustrated in Figs. ^,h. 

As it was noted before, the influence of Qa and flk on the dynamics of Universe 
at early stages was diminished by expansion. This makes the variation of Ac value 
within two arbitrary models smaller (see Fig. Hlb). While we are moving to the 
earlier stages, the values 6c and Ac would converge to the their canonical values 
3/5(37r/2)2/3 ~ 1.686 and IStt^ ~ 178 correspondingly (Fig.|li,b). 

The analytical approximation for the above mentioned value has the following 
form 

Ac = A-B\n{am + C), 

where 

A = 110.89 - 12.91fiA + 7.92{Qa)^ - 0.03/fi^, 
B = -113.28 + 6.84f]A + 7.13(fiA)^ 
C = 0.787 - O.llS^A - 0.093(1]a)^- 

Numerical values of the coefficients are determined using the method of best-fit 
adjustment. Accuracy of the approximation for different models is shown in Fig.Eb- 

5 Mass function of rich clusters of galaxies 

Rich clusters of galaxies are supposed to have been formed at the highest peaks 
of the matter density fluctuations - which, as we know, are almost spherically 
symmetrical. That's why these formations can be described using the spherical 
collapse model. One of the most important features of the large-scale structure 
of the Universe at the scales ~ 10 Mpc is a mass function that determines the 
number of gravitationally bound systems with their masses exceeding the given 
one. To calculate this distribution in the case of spherically symmetrical collapse 
it is convenient to use the Press-Schechter formalism. The amplitude of initial 
perturbations (5 <^ 1) is assumed to be distributed according to the Gauss law 

1 ..^^ 
P(R,t) = ^= e ^^^u^. 
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where 



a\R,T) =< 51 >=^l Sk%\^\W{kR)f 
A fraction of collapsed perturbations equals 

d5 



F(R, t) = / , , -e 



The number of the objects collapsed per unit volume with their masses within dM, 
according to the Press-Schechter approach constitutes ^ 



p dF dR [2 p 4 da{R, 



n(M,r)dM = -2^ — — dM = -A , " , e ^^^(^dM. (28) 

^ ' ^ MdRdM \TTMa^{R,T) dM ^ ' 

The number of objects with masses greater than a given M equals 

N{> M) = / n{M')dM'. 



M 



Dependence of the mass function on cosmological parameters is shown in Fig.[ZK,b. 
For numerical calculations, we have used an analytic form of power spectrum |22| 
normalized according to [2H]. Filled circles depict experimental data taken from 
[28] . Long dashed lines depict the mass functions which differ from the previous ones 
by using canonical value 5c = 1.686 instead of their exact values 5c = 5c{flm,^A)- 
As we can see in Fig. E^-c the difference in the logarithmic scale is small and lesser 
than experimental errors. While switching to linear scales they appear to be more 
prominent. Fig. [HI represents graphs for relative errors of the mass function for 
different cosmological models with fixed values of (a), and for models with zero 
curvature (b) 

N{> M, 5c{nm, Qa)) - N{> M, 5c = 1.686) ^ 
Error = — ^ ^ — — — - — ^ x 100%. 

N{> M,5c{a^,nA)) 

As we can see, the relative difference of the concentrations of rich clusters of galaxies, 
calculated according to the Press-Schechter formalism with real and canonical values 
5c, is quite large and lie within 10 -=-20%. For the closest zero-curvature model with 
Qm = 0.36 and = 0.64 (see Fig. [Zt"), the relative error changes in the range of 
5 to 20 percents. This model is also superior compared to other tests [231 [in]. For 
some other models, this error can be no less than several dozen percents (Fig. [H]). 
In fact, as can be deduced from 



An 

n 



1-4 



A5r 



where A5c is difference between model 5c and its canonical value 1.686. For a{M) < 
0.8 and \A5c/5c\ ~ 0.06 An/n (and correspondingly Error) will be > 20%. This 
fact suggests that one should take into account dependence of critical amplitude 
on cosmological parameters, in order to keep calculation errors within 10%. The 
observational mass function data do not prevent one from possible realization of 
the model with a curvature other than zero. For example, Fig. [HI demonstrates the 
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most consistent with experimental data mass function for the model with a negative 
curvature Qk = 0.1. Theoretical mass functions in the case of elliptical collapse 
Sheth-Tormen (ST) [2E1 and Lee-Shandarin (LS) fI7\ are presented here also. As 
we can see, the theory of spherical collapse well agrees with experimental data, 
since rich clusters of galaxies have formed in maxima of the density peaks which 
predominantly have spherically symmetrical form. 

6 CONCLUSIONS 

We have analyzed the evolution of a spherical dust-like cloud from the early linear 
stage, when the relative difference between the background and densities of cloud 
is (5 -C 1, up to the late non-linear stage, when the collapse occurs, 5 cxd, in 
cosmological models with an arbitrary 3D space curvature and positive values of 
cosmological constant. There are two parameters convenient for comparison of such 
dynamics in various cosmological models. First, the threshold value of initial am- 
plitude of density perturbation, 6c, calculated at the moment of collapse according 
to the law of linear growth. The second is the final amplitude of overdensity fol- 
lowing the dynamical equilibrium, Ac. 6c and are used for calculation of the 
mass function and X-ray temperature of rich clusters of galaxies respectively. On 
the base of the analytical and numerical results presented here it can be concluded 
that for models with negative or zero 3D space curvature and positive cosmological 
constant (A > 0), 6c is smaller than its canonical value 1.686 (fc = A = 0) due to 
longer duration of perturbation increasing. It decreases with increase of flm and 
with decrease of Qk- 

In models with positive 3D curvature 6c can be either greater or smaller in 
comparison with canonical value. In this case critical amplitude 6c goes to infinity 
when flm converges to some small critical value. It happens because the moment of 
slow down of expansion of the background takes place. At this moment amplitude 
of density perturbations grows very quickly tending to the exponential law in models 
with collapse. Such models are gravitationally unstable in the sense that even the 
smallest perturbations cause their collapse in finite time. 

The range of change for the threshold value 6c, within the range of parameters 
-0.4 < fifc < 0.4, < < 1, 0.1 < fi„ < 1 is small: 1.55 < 4 < 1.75. But 
the discrepancy in values of the concentration of rich clusters of galaxies, as calcu- 
lated according to the Press-Schechter formalism, with real and canonical values 6c, 
reaches few dozens of percents for some models. This discrepancy grows with an 
increase of the mass of clusters. 

The range of variation for Ac is much wider, 60 < Ac < 180 (for the same range 
of cosmological parameters Qk, and flm)- This leads to the deviations of the 
calculated values of the X-ray temperature of the gas by ~ 40%, as compared to 
those calculated for the canonical value of Ac, 178. 

Experimental uncertainties of determining mass function and X-ray temperature 
function of rich clusters of galaxies in the vicinity of Milky Way are considerably 
larger than those for the models with Qm ~ 0.4, f^A ~ 0.6 and ~ (see [24] 
and references therein). Nevertheless, with accumulation of data using deep-sky 
surveys (including SDSS, 2dF GRS etc) and with respective increase of accuracy of 
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the data observed, the given dependences of 5c and Ac on VLm and should be 
taken into account when comparing theoretical predictions with the observational 
data in problems aimed searching for adequate models of the Universe. 
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Figure 1: Dependence of the threshold value of the amplitude of density fluctuation 
collapsed at the current epoch ((^c) and calculated for the moment of collapse (a, d), and 
for the moment of turn around of expansion of the spherical cloud (b, e) on the matter 
abundance Q,m for different and f^A. Also, dependence of the age of the Universe for 
corresponding cosmological models expressed in units of the age of the Universe for the 
canonical model with Q,m = 1 and $7^ = VL\ = are presented below (c, f). 
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Figure 2: Dependence of the smallest expansion rate of the Universe (a) and an growth 
factor D{a = 1;^^^) (b) on cosmological parameters. 
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Figure 3: Dependence of the 6c on the scale factor for set of models with zero curvature 
(a), for models with $7^ = —0.4 (b) and $7^ = 0.4 (c). Crosses denote moments of turn 
around of expansion of the perturbation region. 
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Figure 4: Dependence of threshold hnear density fluctuation 6c (a) and overdensity after 
virialization Ac (b) on the moment of collapse Zcoi for different values of the matter density 
and cosmological constant. 
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Figure 5: Dependence of the overdensity after virialization Ac on the matter abundance 
ilm for different fixed values of a cosmological constant (a) and curvature (b). 
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Figure 6: Dependence of the relative error of the analytical approximation for the values 
of 6c (a) and Ag (b) on the parameters of the cosmological model. 
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Figure 7: Mass function of rich clusters of galaxies for models with a zero curvature (a), 
fixed value J7a (b) and the better model with = —0.1 (c) calculated according to the 
Press-Schechter (PS) [T], Sheth-Tormen (ST) (SEI and Lee-Shandarin (LS) ^ h = 0.65 
approaches. 



24 




20 - 





14.6 14.8 15.0 15.2 15.4 
log,„[M/(h-'Mo)] 

Figure 8: Relative difference between concentrations of rich clusters of galaxies calculated 
according to the Press-Schechter formalism. The difference arises from calculations with 
real and canonical values of 6c for fixed values of (a) and for models with a zero- 
curvature (b). 
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